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Abstract. Timely monitoring of the condition of the equipment is an urgent 
task. The paper presents the results of research on methods of monitoring of 
technological equipment, namely, the developed method of determining the in- 
dex of "survivability", which allows evaluating and predicting the condition of 
the unit between maintenance operations. The paper also presents the results of 
the study of monitoring and control systems of technological equipment, name- 
ly, the assessment of the "survivability" indicator and its prediction in the 
STATISTICA program by the exponential smoothing method. 
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1 Introduction 


Large-scale distributed technological objects consist of elements of technological 
equipment combined into technological complexes. 

The absence of a system that includes end-to-end monitoring increases the eco- 
nomic risks of the enterprise, as well as the risks of security breaches, including envi- 
ronmental and energy. Its presence, in turn, allows you to additionally automate the 
ongoing processes, to carry out timely control over changes in performance indica- 
tors. 

Within the framework of this article, the improvement of the system for monitoring 
the technical state of distributed technological objects of oil refineries will be consid- 
ered on the basis of a comprehensive analysis of vibration indicators and the devel- 
opment of a block diagram of the unit for analyzing equipment parameters. 


2 Assessment of the state of technological equipment 


Currently, to assess the reliability of complex technical processes and the types of 
equipment used in them, a large number of working documents, standards [1-2], tech- 
nical conditions and other regulatory frameworks have been circulated. Let us study 
the procedure for carrying out work to assess reliability indicators in accordance with 
RD 50-690-89 [3]: 

a) selection of a reliability test plan; 


b) test planning; 

c) collection of the necessary information; 

d) statistical processing of information. 

The final stage is statistical processing of information to assess the reliability indi- 
cator, carried out by two methods: 

- a nonparametric method with an unknown distribution law, including a direct as- 
sessment of reliability indicators; 

- a parametric method in which the distribution law is known, which also includes 
an estimate of the distribution law parameters included in the formula for calculating 
the estimated reliability indicator, and an estimate of the reliability indicator based on 
the calculated estimates of the distribution law parameters. 

Parametric methods for assessing reliability indicators established in RD 50-690- 
89 are used for exponential, normal, lognormal and diffusion distributions and 
Weibull distribution in the situations indicated in Table 1. 

In the nonparametric method, the problem of determining the distribution function 
and its parameters is not posed. It is necessary to evaluate the values of the reliability 
indicators of the evaluated equipment according to the operational data. 

The disadvantage of the presented technique is its statistical nature. The indicators 
to be assessed are averaged over a large time period, therefore, the forecast with such 
data without taking into account causal relationships and functional analysis of the 
reliability models of objects is approximate 


Table 1. Typical situations for a parametric method for evaluating reliability indica- 














tors 
Designation of 
typical situa- Brief description of typical situations 
tions 
TS-1 Products of the same type are in the same conditions and the 
reliability indicators of these elements are the same 
TS-2 Products of the same type are in different conditions. The ratios 


of the parameters of the distribution laws of elements are known 





Products of the same type are in different conditions. The ranges 
TS-3 of variation of the ratio of the parameters are known depending 
on the location in certain conditions 





Products of the same type are in different conditions. It is a 
TS-4 known fact that the indicators in some conditions are higher 
than in others 














The problem of the lack of representative samples of failures in some cases makes 
it impractical to implement the recommendations of RD 50-690-89 and similar docu- 
ments. 


3 The indicator of ''survivability'' of technological equipment 


The indicator of "survivability" is a dimensionless value that reflects the ability of 
the equipment to work before maintenance while maintaining the trends of the main 
technical characteristics without the occurrence of a certain malfunction [4]. 

The indicator of "survivability" is determined individually for each type of failure, 
since it directly depends on the number of diagnostic parameters, according to which 
the presence of a malfunction is recorded, or rather, on the deviations of real values 
from the target trends. Mathematically, it can be described as follows: 


ae a 
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where w is the deviation of the real value of the parameter from the target; 

i - diagnostic parameter; 

j- 1s the type of failure. 

The following expression is proposed for determining the indicator of 
"survivability" of equipment at time t: 
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4 Determination of the 'survivability'' index of technological 
equipment 


Table 2 presents the diagnostic parameters and their corresponding types of trends 


Table 2. Parameters of the third additional study 























The diagnostic parameter Kind of trend 
Pressure normal 
Bearing temperature escalating 
Vibration velocity of the node normal 








Figure 1 shows a diagram of the change in the "survivability" indicator for the 
study 

Knowing that at time T there is scheduled maintenance, determine the relative 
change in the index "survivability" at time T from its value at time to: 
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Fig. 1. Diagram of changes in the "survivability" indicator for the study 
We make calculations for all the studies and enter the results in Table 3. 


Table 3. Determination of the relative change in the "survivability" parameter 























Study number A, % Ania, % 
1 95,35 
2 93,32 
3 95,48 a8 
4 96,20 











Based on the results of the calculations, we can conclude about the value of the up- 
per limit of the indicator "survivability": when the Ps decreases below 95% of the 
maximum value it is necessary to carry out minimum maintenance of equipment (in- 
spection, check leaks, the appearance of noise and other characteristic markers of 
failure). 

To determine the lower boundary of the Ps, at reaching which an urgent unsched- 
uled repair is necessary, it is necessary to use the prediction of the resulting time se- 
ries. 


5 Predicting the ''survivability" index 


We select several key parameters, with the help of which we will further predict 
the probability of failure until the nearest maintenance: 

- vibration characteristics (in our case, vibration velocity); 

- the head produced by the pump; 

- bearing temperature. 

We will use the data samples obtained from the database of artificially simulated 
time series to construct graphs of indicators change in time [5]. 





The production head of the pump varies in the interval [50-5%; 50+7%]; the max- 
imum bearing temperature must not exceed the room temperature by more than 50°C 
and the set limit value of 363 K (90°C); the limit value of the vibration velocity of 
some units of the equipment does not exceed 11.2 mm/sec and cannot exceed the 
value of the indicator equal to 7.1 mm/sec for more than 30 days. 

Samples consist of 60 indicators of a certain parameter, recorded every 12 hours. 
The choice of the time range is due to the manufacturer's recommended maintenance 
period [6], and is about 700-730 hours. 

As a result of construction of diagrams of changes of parameters in time, as well as 
limit and target trend lines [7], the following graphs were obtained (Fig. 2, 3, 4). 

After calculating the total maximum deviations, we determine the value of the de- 
viation at each moment of time, summarize these deviations, and then normalize each 
of the indicators by the following formula: 


Ww, = ees, (3) 
(>) AO), 


where w; — is the normalized deviation; 
Aq, - the value of the deviation at time f; 


ey AQ), - is the value of the total maximum deviation from the trend at time ¢. 


We calculate the first and second terms of equation 3, given the values of the tun- 
ing coefficients k, m, and 4. After that, we determine the survivability index of the 
equipment for the initially proposed type of failure. 

The change in the survivability index over time is shown in Figure 5. 

In the period from 0 to 84 hours, the survivability index decreased, which may in- 
dicate the equipment running-in phase. Any deviations of diagnostic parameters dur- 
ing this time show that the equipment enters normal operation mode, and its parame- 
ters are only approaching the trend lines. 
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Fig. 2. Diagram of head changes of pump K 80 


— Temperature value °c 


—— Temperature trend line 


Ho Temperature limit 









































































































































Fig. 3. Temperature variation diagram of the K 80-50-200-E pump bearing 
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Fig. 4. Vibration velocity diagram of the investigated pump unit K 80 
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Fig. 5. Diagram of change in the "survivability" index of a centrifugal pump for one 
type of failures 


Then the value starts to rise again and takes on values slightly deviating from the 
average value of approximately 0.950. This level is maintained for the remainder of 
the life of the pump before maintenance. 

On the basis of one-time conducted research it is impossible to judge precisely 
about boundary values of "survivability" parameter, at which minimum maintenance 
up to planned maintenance is required, or equipment is taken out of operation for 
urgent repair, restoration of damaged unit and bringing of diagnostic parameters to the 
trend values. Therefore, it is proposed to conduct at least three additional studies, 
using the same diagnostic parameters, but changing the data samples. At the same 
time, the target trends of indicators of head, bearing temperature and vibration veloci- 
ty of the unit remain the same, since we will assume that each of the studies was car- 
ried out on the same type of equipment 

Next, we will use the software package STATISTICA to predict the index of "sur- 
vivability”. 

Short-term or medium-term forecasting will be preferable in the conditions of the 
analysis unit. The most suitable forecasting methods are exponential smoothing and 
autoregressive integrated moving average (ARIMA). 

Let us perform the forecasting of RL with the help of exponential smoothing. 

To build the prediction, we choose the change in the "survivability" index obtained 
from the results of the first additional study, because the diagram (Fig. 6) shows the 
decrease in RR over time. 

In STATISTICA program, we turn to the Time Series/Forecasting module. The 
next step is to assign the variable under study (Variables) - the index of "survivorship" 
from the original study. 

Next, you should select the Exponential smoothing & forecasting item. This item 
offers 12 different variants of exponential smoothing models. For our case, the expo- 
nential model will be suitable. 


To perform exponential smoothing for the selected model it is required to define 
two smoothing parameters - Alpha (a) and Gamma (7). 

For the smoothing parameter a, which is responsible for smoothing of the observa- 
tions directly, the following is true: if a is 1, then the previous observations are com- 
pletely ignored; if a is 0, then the current observations are ignored [8-9]. 

For the smoothing parameter y, which is responsible for smoothing of the trend: If 
y equals 0, the trend is constant for all values of the time series (and for all forecasts); 
if y equals 1, the trend is determined mostly by observational errors. 

We will search for the parameters using the Grid search tab. Here we set start, step 
and limit values for the parameters a and y. Define these values, also select the item 
Display parameters for 10 smallest mean squares and click Perform grid search. 

The window that opens displays a table with the 10 best models. The determining 
factor for model selection is the Mean Absolute Error shown in the outermost column. 


Parameter grid search (Smallest abs. errors are highlighted) (NporHosuposanue.sta) 
Model: Expon. trend, no season ; S0=1,012 TO0=,9702 
VUccnegosanne Not 


Model Alpha feces Mean | Mean Abs | Sums of | Mean | Mean % |Mean Abs 
Number | Error Error | Squares | Squares | Error % Error 
1850 0.72100010,981000 0.000677 0.001771 0.000682 0.000011 0.070104 0,186003 









1900 0.741000 0,981000 0,000659 0,001767 0.000683 0,000011 0.068220 0,185543 
1800 0,701000 0.981000 0,000696 0,001778 0,000684 0,000011 0.072091 0,186834 
1899 0.741000 0.961000 0,000673 0,001771 0,000686 0.000011 0.069642 0,185968 
1849 0,721000 0.961000 0.000691 0,001776 0,000686 0.000011 0.071565 0,186524 
1950 0.761000 0.981000 0.000642 0001765 0,000687 0,000011 0.066429 0.185362 
1750 0,681000 0.981000 0,000716 0,001790 0,000689 0,000011 0.074189 0,188112 
1949 0.761000 0.961000 0,000655 0,001768 0,000689 0.000011 0.067815 0.185688 
1799 0.701000 0.961000 0.000711 0,001785 0,000689 0.000011 0.073593 0.187501 
1898 0.741000 0.941000 0,000687 0,001775 0,000690 0.000011 0.071126 0.186440 


Fig. 6. Search result on the grid 


The model with the lowest mean absolute error (1950) has the following parame- 


ters o and y: 
a= 0,761; 
y = 0,981. 


We select 14 periods for the forecast, which will be equal to 168 hours of equip- 
ment operation or one week. Then click OK (Summary) and we get the graph shown 


in Figure 7. 
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The blue color on the graph shows the real values of the "survivability" indicator, 
the red dotted line - the forecast that is obtained by using exponential smoothing. 


The green dotted line is the residuals corr 
actual values and the values obtained by smo 
of operation of the equipment do not exceed 4 

It follows from the graph that if the equip 





esponding to the difference between the 
othing. The residuals for a longer period 
+0.0025. 

ment continues to operate after the time 


of planned maintenance without its implementation, the indicator "survivability" will 
continue to decline and after 168 hours will reach the value of 0.927 (Table 4). 


Table 4. Prediction results with exponential s 


moothing 





















































Period bee oe aie The value of the "survivability" indicator 
61 720 0,929828 
62 732 0,929680 
63 744 0,929532 
64 756 0,929384 
65 768 0,929236 
66 780 0,929089 
67 792 0,928941 
68 804 0,928793 
69 816 0,928645 
70 828 0,928498 
71 840 0,928350 
72 852 0,928202 
73 864 0,928055 
74 876 0,927907 
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In order to compare the results, we make a prediction using the ARIMA method. 

In the STATISTICA program we access the module Time Series/Forecasting. The 
next step is to select the variable under study (Variables) - the index of "survivability" 
from the initial study and construct a linear graph for it. Press the OK button and in 
the dialog box that opens select the Plot of the selected variable in the Review & plot 
tab. 

For further verification of the time series it is necessary to perform the autocorrela- 
tion analysis. 

In the previous dialog box, select the Autocorrelations tab (Autocorrs). We set the 
number of lags to 14, which in our case will correspond to 14 periods of measure- 
ments or 168 hours of equipment operation (one week). Then we build an autocorrela- 
tion function (Autocorrelations). 

The resulting autocorrelation function is shown in Figure 8. 


Autocorrelation Function 


Study #1 


(Standard errors are white-noise estimates) 














0 ~~~ Conf. Limit 





aoe | 


-1,0 05 


o 


Fig. 8. Autocorrelation function 


According to the results of the autocorrelation function, the first lag (measurement 
period) is highly significant and significantly exceeds the limits of the confidence 
interval. 

Let us remove the significant value of autocorrelation by taking the difference of 
the order of 1. Return to the analysis dialog box and in the Difference, Integrate tab 
assign the value of lags equal to | in the Difference item. After that press OK and we 
obtain the transformed time series. To check if the autocorrelation value was changed 
at the first lag, it is necessary to plot the autocorrelation function for the transformed 
time series (Figure 9). 
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Autocorrelation Function 
Study #1; D(-1) 
(Standard errors are white-noise estimates) 












Lag Corr. S.E. Q p 
1 +,300 ,1269} , (ae) 5,60 ,0180 
2 +,243 ,1258 } : (saa 9,33 ,0094 
3 +,125 ,1247} i Bi 10,33 ,0159 
4  +,070 ,1236} {_g- 10,66 ,0307 
5 +,128 ,1225} : £8: 11,75 ,0385 
6 -,030 ,1214} i ] 11,81 ,0665 
7 -,043 ,1202} : I} : 11,93 ,1029 
8 -,060 ,1190} 0 12,18 ,1433 
9 -,095 ,1179} : | : 12,83 ,1707 
10 -,065 ,1167} : 0 : 13,14 ,2160 
11 -,162 ,1155} | : 15,10 ,1780 
12 -,082 ,1143} : a 15,62 ,2093 
13. -,044 ,1131} ' H} ' 15,77 ,2617 
14 -,029 ,1118} : ] : 15,84 ,3234 
ow — —— 0 ~ Conf. Limit 
-1,0 -0,5 0,0 0.5 1,0 


Fig. 9. Autocorrelation function of the transformed time series 


Taking the difference with lag 1 led to a decrease in the autocorrelation function on 
the corresponding lag. Also, based on the obtained graph, we confirm the absence of 
any periodicity (seasonality). 

To select the parameters of the ARPSS model, it is necessary to construct two 
functions, the autocorrelation function and the partial autocorrelation function for the 
transformed series. 

To build the partial autocorrelation function, select the item Partial Autocorrela- 
tions in the tab Autocorrelations. 

Figure 10 shows the resulting partial autocorrelation function. 


Partial Autocorrelation Function 
Study #1: D(-1) 
(Standard errors assume AR order of k-1) 
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Fig. 10. Partial autocorrelation function of the transformed time series 
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For the partial autocorrelation we observe the following: no significant lags, no 
seasonality, the graph decays exponentially. 

According to the rules for determining the ARPSS parameters, we draw the follow- 
ing conclusions: 

- if the autocorrelation function decreases exponentially and the partial autocorrela- 
tion function has a sharply prominent value for lag 1, and there are no correlations at 
other lags, then the p value corresponding to autoregression (Autoregression, AR) is 
1; 

- If the autocorrelation function has a sharply distinctive value at lag 1, there are no 
correlations at other lags, and the partial autocorrelation function exponentially de- 
creases, then the index q corresponding to the Moving Average (MA) equals 1. 

Let's return to the time series analysis menu. In order to estimate the prediction er- 
ror, we select a part of the data from the sample for cross-checking. To do this, it will 
be necessary to exclude a certain range of data from the time series. In the Condition 
item we set the range of observations from 1 to 46. Then 47-60 observations will 
remain for cross-checking. 

In the time series analysis menu, select ARIMA & autocorrelation functions. In the 
dialog box that opens, set the parameters that were previously defined. 

In the ARIMA-model menu, select the Advanced tab. In order to obtain numerical 
values of the predicted data, as well as the graph of the forecast, let's set: 

- number of observations (number of cases) - 14; 

- start at case - 47; 

- confidence interval - 0,65. 

We decrease the confidence interval from 90% to 65% in order to narrow the fore- 
casting range, which will make it possible to obtain more accurate forecast values. 

Then select the items Forecast and Plot series & forecasts. The numerical results of 
the forecast and their graphical representation are shown below (Figure 11). 


Forecasts; Model:(1,1,1) Seasonal lag: 12 
Input: Study #1 
Start of origin: 1 End of origin: 46 
Include cases: 1:46 
1,02 : - ° - - - = - ~ . . . . 1,02 





1,00 1,00 
0,98 0,98 
0,96 0,96 


0,94 0,94 











0,92 0,92 
0,90 0,90 
0,88 0,88 
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— Observed -- Forecast + 65,0000% 


Fig. 11. The graph of the "survivability" indicator with the ARIMA model verifica- 
tion 
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Let's perform cross-check to determine the value of absolute error of predicted val- 
ues of time series. The results of cross-check are presented in Table 5. 


Table 5. The result of cross-checking the predicted series 





















































No. Observable Beeaiciea tories The magnitude of the absolute 
series error 
47 0,9311 0,931599 0,00047 
48 0,9313 0,931727 0,00035 
49 0,9303 0,931816 0,0015 
50 0,9297 0,931879 0,0023 
51 0,9296 0,931922 0,0024 
52 0,9291 0,931953 0,0030 
53 0,9289 0,931975 0,0032 
54 0,9296 0,931990 0,0025 
55 0,9298 0,932001 0,0023 
56 0,9305 0,932008 0,0016 
57 0,9298 0,932013 0,0023 
58 0,9302 0,932017 0,0019 
59 0,9299 0,932020 0,0022 
60 0,9311 0,931599 0,00047 














The magnitude of the error is extremely small, which allows us to conclude that 
this prediction method is also suitable for use in the parameter analysis block. 


6 Conclusion 


The study examined one of the key tasks of rational operation of refinery units - 

monitoring the condition of equipment. Untimely monitoring of increasing operation- 
al and functional indicators of units leads to sudden breakdowns and accidents. 
We developed an analysis unit based on determining the parameter "survivorship" - a 
dimensionless value, reflecting the ability of the equipment to keep working until the 
moment of maintenance, as well as predict the possibility of operating the unit with- 
out maintenance, while maintaining the same level of trends in the main diagnostic 
parameters. 

Use of complex analysis of parameters, characterizing state of equipment, will al- 
low personnel to receive timely signals about occurrence of failures and to eliminate 
them, avoiding damage of technological equipment, disturbance of technological 
process and providing safety of personnel. These measures, in case of sudden break- 
age or accident, make it possible to prevent economic losses of the enterprise. 

Due to timely control, it becomes possible to predict the condition of the equip- 
ment before and after the scheduled maintenance. For this purpose, the method of 
forecasting the index of "survivability" of technological equipment in the 
STATISTICA program was considered and developed. 
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